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Yang-Mills classical solutions and fermionic zero modes from lattice 
calculations * 
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We study a series of problems in classical Yang-Mills theories using lattice methods. We first investigate SU(N) 
self-dual configurations on the torus with twisted boundary conditions. We also study the zero modes of the Dirac 
equation in topological non-trivial background Yang-Mills fields. 



1. Introduction 

It has long been argued that classical configu- 
rations might play a relevant role in the under- 
standing of low-energy phenomena like confine- 
ment and chiral symmetry breaking. Here we will 
present a study of a particular kind of solutions 
which appear on the torus when twisted bound- 
ary conditions are imposed. These are instanton- 
likc configurations with the peculiarity of having 
fractional topological charge. In the large N limit 
their action scales as 8-k 2 /Ng 2 . As has been ar- 
gued in 0| they may be important to understand 
the long-distance dynamics of QCD. In this pa- 
per we will present a preliminary investigation 
of: i) the large N limit behaviour ii) the zero 
modes of the Dirac operator on the background 
of fractional-charge instantons. 

2. Yang-Mills classical solutions 

The first part of this work is a study of SU (N) 
(anti) self-dual Yang-Mills classical solutions on 
the torus with twisted boundary conditions. It is 
known that in this case the topological charge Q 
is not always an integer H, but is given by 



2 / Tr [ F ^ F ^)d i x = v-—, (1) 
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and F^y its dual. Here 



with n^ v the usual twist tensor and ki 
n i: j = e ijk m k . 

From Schwarz-inequality 

S =lJ Tr ( F ^ F ^) d±x> 8n 2 \Q\ 



(2) 
n 0i , 

(3) 



one easily derives that if n ^ (modulo N ) an ob- 
struction is found for zero-action configurations. 
We are interested in those solutions with minimal 
non-trivial action 



8f_ 
N ' 



(4) 



where v and K are integers, F^y is the Yang-Mills 
strength tensor in the fundamental representation 



on a volume [0, L} 3 x [0,T], with T > L. Some 
of them are already known: i) For specific values 
of N, 't Hooft has constructed non-abelian so- 
lutions with constant field strength which turn 
out to be (anti)self-dual whenever the sides of 
the torus satisfy certain relations (see || for de- 
tails), ii) there are also numerical studies of the 
solution for SU(2) with twist m = (1,1,1) and 
k = (1,1,1) 0. The present work is inspired by 
this last reference. 

We have restricted our analysis to the following 
twist tensors: 

1. Spatial twist, always fh = (1, 1, 1). 
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2. Temporal twist, two cases: 
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• k = (1,0,0) for N = 3,4, ,...,13 the 
solution is in this case antiself-dual, 
Q = -l/N. 

• k = (n, n, n) for N = 3n + 1 = 4, 7, 
10, 13 the solution is here self-dual, 
Q = l/N. 

The solutions were generated on Nf x N t lat- 
tices (jV t >> N s ) by cooling with the Cabibbo- 
Marinari-Okawa algorithm [||. They verify the 
following properties: 

1. Scaling towards the continuum solution. 

2. Self-duality or antiself-duality. 

3. The energy profile 

e(t) = j S(x,t)d 3 x (5) 

with S(x, t) the action density, is instan- 
tonic and independent of temporal twist. 

4. Behavior for large N: 

• The energy profile behaves as e(t) ~ 
f(t/N)/N 2 as illustrated in Figure 1. 

• The norm of the field strength tensor 
becomes independent of the spatial co- 
ordinates. 

• In the gauge 

A = A t {t = -oo,x) = 0, (6) 

one eigenvalue (the same one for A4, 
Bi, Ei) dominates the solution in the 
sense that it gives the major contribu- 
tion to the action. 

A more detailed description of the properties 
of these solutions will be presented in |(| . 

3. Fermionic zero modes 

The second part of this work is a study of 
fermionic zero modes of the Dirac equation in 
the background of topologically non-trivial Yang- 
Mills fields. 

To obtain the zero modes we have used the lat- 
tice approach with Wilson fcrmions and also with 




Figure 1. The energy profile for 4 different 
SU(N) solutions multiplied by iV 2 is plotted as a 
function of t/N. 



naive (Kogut-Susskind) fcrmions. In the massless 
case, the Dirac operator reads 

W«) = -\ X)K r - i»)M n Mn + A) + 

(r + 7^)U^(n - fi)^{n - /})] + 4r*(n), 

where \t is the fermion field, the gauge field, 
7 M the Dirac matrices and r the Wilson parame- 
ter. 

We look for the lowest eigenvalues (and the cor- 
responding eigenvectors) of the operator (75 p) 2 
using three different methods: 

• Conjugate gradient 

• Local minimization 

• Lanczos (only for eigenvalues) 

The results we have obtained are the same inde- 
pendently of the method employed in the diago- 
nalization. We should point out that the Wilson 
parameter must be big enough to avoid fermionic 
doubling and, at the same time, small enough 
such as not to distort the solution. For a lattice 
size with N s > 8 a good value is r = 0.001. 
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Table 1 

The three lowest eigenvalues of (75 ft) 2 in the 
presence of a Q = 1/2 instanton on a Nf lattice. 



N s 


l s *(xlCT 9 ) 


2™ d (xicr y ) 


3 rd (xi(r 6 ) 


11 


6.946 


6.946 


4.150 


12 


4.835 


4.835 


4.140 


13 


3.546 


3.546 


4.107 


14 


2.633 


2.633 


4.103 



As a check of the method we successfully ex- 
tracted the zero mode for the instanton solution, 
also known as the 't Hooft zero mode. 

Here, we will present our results for the zero 
modes in the background of the SU(2) self-dual 
gauge field with Q = 1/2 on the [0, L] 4 torus with 
twist m = k = (1,1,1). Fermions are taken in the 
adjoint representation of SU(2) (to avoid twist 
singularities). 

We have worked on Nf lattices with N s = 
11,12,13,14 and r = 0.001. The three lowest 
eigenvalues obtained are shown in Table 1. We see 
that two of these eigenvalues scale to zero. These 
would correspond to the two zero modes predicted 
by the Index Theorem. They are related to each 
other by the operation $2 = 7173^1 • Our results 
show that the gauge invariant density of these 
lowest eigenvectors scales nicely towards a con- 
tinuum curve. Furthermore, this density is cu- 
bic symmetric and peaks at the position of the 
Q = 1/2 instanton. This is examplified by Figure 
2, where we display the quantity: 



*(ar 4 ) 



E 



^ai*^ai dxidx2dx3 



(7) 



In the above equation, £4 is one of {x, y, z, f}, and 
Xi,X2 ) %3 are the other three components, where 
the origin of coordinates is taken at the maximum 
of the gauge field configuration. It is clear from 
the figure the nice scaling of the results. The 
curves for naive fermions are not displayed, but 
fall on top of the other. 

A more detailed description of these solutions 
and of others obtained for different classical Yang- 
Mills backgrounds will be given in M. 
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Figure 2. The quantity $ is plotted as a function 
of X4 for different lattice sizes. 



for MGP is thankfully acknowledged. AGA and 
AM would like to thank CICYT (Grant AEN 
96 — 1664) for their financial support. 

REFERENCES 

1. A. Gonzalez- Arroyo and A. Montero, Phys. 
Lett. B387 (1996), 823. 

2. G. 't Hooft, Nucl. Phys. B153 (1979), 141. G. 
't Hooft, Acta Phys. Austriaca, Suppl. XXII 
(1980), 531. 

3. G. 't Hooft, Commun. Math. Phys. 81 (1981), 
267. 

4. M. Garcia Perez, A. Gonzalez-Arroyo and 
B. Soderberg, Phys. Lett. B235 (1990), 117. 
M. Garcia Perez and A. Gonzalez-Arroyo, J. 
Phys. A26 (1993), 2667. 

5. N. Cabibbo and E. Marinari, Phys. Lett. 
B119, (1982), 387. M. Okawa, Phys. Rev. 
Lett. 49 (1982), 353. 

6. M. Garcia Perez, A. Gonzalez-Arroyo and A. 
Montero, in preparation. 

7. A. Gonzalez-Arroyo, A. Montero and C. 
Pena, in preparation. 



Acknowledgments: DFG financial support 



